Nonequilibrium propagation and decay of a bound pair in driven t—J models 
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We perform an accurate time-dependent numerical study of out-of-equilibrium response of a 
bound state within t-J systems on a two-leg ladder and a square lattice. We show that the bound 
hole pair decays with the onset of finite steady current if both mechanisms for binding and the dis- 
sipation share matching degrees of freedom. Moreover, by investigating the mechanism of decay on 
the square lattice we find that the dynamics is governed by the decay in the direction perpendicular 
to the electric field, leading to much shorter decay times in comparison to the ladder where such 
dynamics is topologically restricted. 

PACS numbers: 71.27.+a, 71.38.Mx, 74.20.Mn, 74.40.Gh 



Introduction. — Equilibrium states of correlated many- 
body systems can exhibit many distinct collective phe- 
nomena, with the Mott insulating, superconducting (SC) 
and magnetically ordered states being the most promi- 
nent examples. A fundamental question in this respect 
is how such many-body states behave when driven far 
from their equilibrium by finite external fields. Many 
theoretical studies have recently focused on the dielectric 
breakdown of the Mott insulator [l|, Q and the majority 
of them addressed this problem within the framework 
of the half-filled Hubbard model. Much less is known 
about nonequilibrium properties of strongly-correlated 
systems away from half-filling, where unconventional su- 
perconductivity may possibly emerge. Ballistic response 
of superconductors to constant electric field has recently 
been observed It shows up as a current that in- 
creases linearly in time up to a threshold value above 
which the superconductivity is destroyed. The properties 
of cuprate superconductors subject to the electric current 
have also been studied within a variational (equilibrium) 
approach For low concentration of holes, sufficiently 
strong current destroys superfluid stiffness while pairing 
remains intact. 

In this Letter, we do not discuss superconductivity. In- 
stead we focus on the real-time nonequilibrium dynam- 
ics of its basic ingredients, i.e., bound pairs of charge 
carriers. According to the Ohm's law, driving of charge 
carriers by constant electric field leads to a finite time- 
independent current due to a steady emission of excita- 
tions, e.g., phonons or magnons. Here, we address a ques- 
tion whether a bound pair of two carriers can respond in 
a similar way, i.e., whether it can acquire a constant ve- 
locity upon constant electric field without decaying into 
two separate carriers. In principle, such problem can be 
addressed in experiments on systems in which pairing 
precedes the SC phase coherence. 



swer is not immediately obvious since pairing degrees 
of freedom simultaneously present quite effective dissi- 
pation channels. Therefore, propagation of a bound pair 
under electric field causes steady emission (heating) of 
magnons or spin fluctuations which simultaneously me- 
diate the pairing interaction. Here, we consider two holes 
in the <-J ladder and square lattice and carry out fully 
microscopic calculations taking into account two most 
relevant phenomena: (i) pairing by the exchange of spin 
excitations that has been so far studied predominantly 
under equilibrium conditions 043 i dissipation by 
emission of spin excitations Q. Recent investigations of 
driven systems at half-filling have contributed to a gen- 
eral understanding of heating in isolated systems which 
on a long time scale suppresses any steady current [ol-fl^. 
There are two possible ways to avoid this problem: ei- 
ther to couple the system to enivironment |lll4l5| . or to 
consider vanishingly small concentration of charge carri- 
ers d, [l^ . In our work, we consider the second option. 
For a fixed number of carriers (two holes) heating is a fi- 
nite size effect Q , thus we restrict our analysis to regimes 
where the results are essentially size independent. 

Another important property of driven quantum sys- 
tems that remains a widely unexplored subject is the 
role of dimensionality transpiring in the compelling phe- 
nomena emerging in the direction perpendicular to the 
driving [l^ [13] • The most common numerical ap- 
proaches to nonequilibrium correlated systems were de- 
veloped for studies of either one-dimensional [isf or 
high-dimensional 3 23 systems. By using dynamical 
mean-field theory it has been shown that in the limit of 
extremely strong electric field, a D-dimensional system 
exhibits equilibrium properties in _D — 1 dimensions per- 
pendicular to the field [13|. We show that for moderate 
electric field an unexpected effect emerges in the 2D sys- 
tem: the decay of the bound state is governed by the 



Despite a seeming simplicity of this problem the an- motion of charge carriers perpendicular to the field. 
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Model and setup. — We consider a driven t-J model 
with two holes (also referred to as charge carriers) on 
the ladder and the square lattice 



-H.c.)+^Jij(SiSj-ininj) 



(ij> 



(1) 

where Ci,s = Ci,s(l — rn^-s) is a projected fermion op- 
erator, to represents nearest neighbor overlap integral, 
the sum (ij) runs over pairs of nearest neighbors and 
hi = ni^t' + '^i,^ ~ '^'^■1,^1^1,1 is a projected electron number 
operator. On the ladder, Jy may be different for interac- 
tions along and perpendicular to the ladder's leg, while 
we set Jjj = J for the square lattice. The constant elec- 
tric field F is switched on at t = 0. It is applied along the 
ladder's leg and along the diagonal of the square lattice, 
i.e., we set = -i+)Ft and = -(-h)Fi/v^ 

for positive (negative) directions of carrier hopping in 
Eq. (HI, respectively [2ll |. 

We apply the time-dependent exact diagonalization 
method (t-ED) within the full Hilbert space to calculate 
out-of-equilibrium response of the driven t-J ladder with 
periodic boundary conditions, while we use the time- 
dependent exact diagonalization method defined over a 
limited functional space (t-EDLFS) for the t-J square 
lattice [1, [2^ . The latter method has been successfully 
applied to calculation of the ground state of the t-J 
model with two doped holes ; the details of the method 
are given elsewhere Q • 

Time evolution of both systems subjected to an exter- 
nal electric field is calculated using the iterative Lanc- 
zos method [2J]. This method has been lately ap- 



plied to calculate out-of equilibrium response of differ- 
ent quantum many-body systems, both during the con- 
stant driving [ig as well as during and after photoex- 
citations [25|. In contrast to the large part of recent 
nonecruilibrium studies of driven strongly-correlated sys- 
tems [2|,H,[l3, 3 ill I ou'^ approach enables calculation of 



the steady state where the driven charge carriers acquire 
constant velocity due to the propagation in dissipative 
medium. The latter, which may also be referred to as a 
quantum heat bath, can be either modeled by interaction 
with magnons phonons 27|, or both Interest- 
ingly, properties of these systems share some similarities 
with driven systems at half-filling, with the most promi- 
nent example being the current-field (j-F) characteris- 



tics 0,1, [Hill [11 



28| . However, the calculation oij-F 
characteristics is not the main goal of the present study. 
Instead, we rather focus on the conditions and mecha- 
nism of decay of the bound state under the influence of 
the electric field. 

Propagation of a hound state on the ladder. — The re- 
sponse of a quantum system to a constant electric field 
considerably depends on the strength of the field. For 
very small F ^ the adiabatic regime (AR) with zero 
net current is observed, while larger F gives rise to 



the dissipative regime (PR) where constant F induces 
a finite dc current 0, Throughout the work, we 

introduce the average distance between holes in par- 
allel ((i||) and perpendicular (d^) direction, defined as 

d^. = v'Er (i--e^.)'C(r) - d^{t = 0), where C(r) = 
^.(n('ri,[Yr)/2 [n^ is a hole number operator), while e|| 
and ej^ are unit vectors parallel and perpendicular to the 
field, respectively. 

The existence of the bound pairs in AR is rather un- 
ambiguous. If the spin gap (between the ground state 
and the excited states) is finite for arbitrary momentum, 
then the system remains in its instantaneous eigcnstate. 
The distance between carriers c?|| {t) and the change of the 
total energy E{t) arc determined by the instantaneous 
value of the adiabatic parameter (t){t) = Ft. Since the 
Hamiltonian H[4)] is periodic, d\\ (t) and E{t) are periodic 
as well. They oscillate with a frequency twice larger than 
the Bloch frequency lob = F, where doubling originates 
from the (double) charge of the bound pair. Plotting the 
system evolution in the plane [E{t),d\\[t)\, closed loops 
emerge as a result of periodicity. Such a loop is shown in 
Fig. [Ija), see the (red) curve for the weakest field. The 
horizontal extent of this loop is determined by the dis- 
persion of the ground state. Here, we show results for 
anisotropic J when exchange interaction along the rungs 
J± is much stronger than the interaction along the legs 
J|| . In such case the spin gap occurs for arbitrary flux 
and the presence of AR for F ^ is unquestionable. 




Figure 1. (Color online) Time evolution on a 2 x 10 ladder 
with two holes under constant F switched on at f = 0. (a) 
d|| it) vs E{t) for Jx = 1 (along the rungs), Jy = 0.4 (along the 
legs) and various F. (b) x{t) (main) and d\\{t) vs x(t) (inset) 
for three different values of V and W . We use J± = Jy = 0.4 
and F = 0.125 in (b). In the inset, d||(t) is not subtracted 
by d||(t = 0); thin dashed horizontal line shows dy for two 
noninteracting fermions. 

Contrary to AR, DR is characterized by a steady in- 
crease of energy E{t). Therefore, moving from AR to 
DR must be accompanied by a destruction of the loops 
in the [E{t),d\\{t)\ plane. The direction of the loop de- 
formation shows whether bound pairs survive in DR. If 
they do, then the loops should be elongated only in the 
horizontal direction [E[t) increases], while in the verti- 
cal direction d\\{t) should remain bound. However, as 
shown in Fig. [IJa) the opposite happens: the increase of 
fi|| {t) goes along with the increase of energy. This indi- 
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cates that the bound pairs dissociate immediately when 
DR sets in. Carrying out calculations for various parame- 
ters and fields, we found no case with propagating bound 
pairs in DR. 

The central question is why is it so and what mecha- 
nism would allow pairs to propagate with a steady veloc- 
ity in a dissipative environment. The propagation with 
a steady velocity under a constant F leads to a steady 
increase of energy. Since the average kinetic energy of 
charge carriers should remain constant (steady current 
in DR), the electrostatic energy has to be transformed 
into excitations of the spin background. However, the 
interaction of holes with the spin background is simulta- 
neously the only pairing mechanism. 

Below we demonstrate that this double role of spin 
excitations - as a dissipation mechanism and as a pair- 
ing mechanism - is responsible for decay of bound pairs 
in DR. For this sake we extend the t-J model by the 
nearest {V) and next nearest neighbor [W) attractive 
interactions 

(ij> 

which play the role of additional pairing mechanisms. 
Fig. [Ijb) demonstrates how this mechanism affects d|| (t) 
and the distance travelled by the center of mass x{t) for 
various V but constant V ~ W. We notice that chang- 
ing of V and W does not influence x{t) and the steady 
increase of energy E{t) = 2Fx{t) being a hallmark of 
DR is clearly visible. While in the pure t~J model the 
pairs dissociate very quickly, see the curve V = W = 
in the inset of Fig. [TJb), both holes stay together for suf- 
ficiently large V and W. Therefore, we notice that the 
bound pair of carriers can propagate under constant F 
with a steady velocity provided that there are different 
mechanisms (degrees of freedom) responsible for pairing 
and dissipation. Introducing simultaneously V and W 
leads to the attractive potential between holes that allows 
hopping of the hole pair without breaking the attractive 
potential. The opposite case of large V and W = leads 
to pair breaking. 

Inclusion of the interactions V and W allowed us to 
explain why bound pairs in the pure t-J model (an pos- 
sibly also in the electron-phonon systems decay as 
soon as the response is dominated by the dc current. We 
would like to emphasize though that replacing attractive 
V and W terms with more realistic electron-phonon in- 
teraction would introduce additional channel for the dis- 
sipation of potential energy through emission of phonons. 
It is thus plausible to speculate, that a spin-lattice bound 
pair would as well decay in the dissipative regime. From 
now on, we again focus on the isotropic t-J model and 
setV = W = 0. 

Mechanism of decay. — We now extend out study to the 
square lattice. We compare the response of the system 
in the parallel vs perpendicular direction with respect to 



the electric field. When investigating properties of the 
bound state, see Fig. [2|a), we observe d±{t) > d\\{t) for 
all times of our calculation. This implies, in contrast to 
the ladder system, that the dynamics perpendicular to 
the electric field governs the decay of the bound state for 
(at least) short and intermediate times. We focus on this 
issue further on to clarify the mechanism of decay on the 
square lattice. 




Figure 2. (Color online) Driven hole pair on the square lattice, 
(a) d||(t) and d±{t) for J = 0.4 and F = 0.7. (b) d±{t)J vs 
x(t),J for different J and F. Black solid line represents a fit 
dx(t) = /3a::(t)exp(-^^j|j^), where P = 1.474 and a = 2.066. 
(c) Current j{t) for different J and F. (d) ^(t) vs x{t), see 
Eq. ([6]), for the same set of parameters as in (b). 

We first investigate the transient time t < t* , where 
t* is defined as the characteristic time needed to reach 
the steady state after turning on the field. We may ex- 
pect that the dynamics of decay at t < t* is strongly de- 
pendent on properties of the bound state in equilibrium, 
which is determined by the energy scale J Q- Remark- 
ably, numerical data reveal both J- and i^-independent 
scaling of d±{t) vs x{t) at short times. As shown in 
Fig. m^b), the universal relation between d±{t) and x{t) 
can be well described by 

dj_{t) (xx{t) e~^^^, (3) 

with a ^ 2. Such relation, which indicates an activated- 
type of behavior, incorporates the information about the 
size of the bound state at t = 0. Indeed, the average dis- 
tance between two holes in the equilibrium scales roughly 
with 1/J for J > 0.4 [111. For i > t* the scaling of 
Eq. breaks down, however, d±{t) vs x{t) still show a 
_F-independent behavior for a fixed J. In the following, 
we will show that such _F-independent scaling is as well 
expected for larger times t* . 

We now proceed to describe the properties of the 
steady state, which is defined as the regime when the cur- 
rent along the field j{t) is time independent. In Fig.[5Jc) 
we plot j{t) which clearly marks the onset of constant 
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Figure 3. (Color online) C(r) measuring time-dependent probability for the hole pair to be at a relative position r, fulfilling 
the sum rule X!j.C(r) = 1. Arrows mark direction of electric field. We set J = 0.4 and F = 0.7. 



current for t > 2. In this regime holes are already well 
separated in space and the binding mechanism should be 
significantly weaker than in the equilibrium. Therefore, 
one may expect a rather independent movement of holes 
in the direction perpendicular to F. This movement can 
be naively modeled by a random-walking process: 



(4) 



where (t~t*)D± is proportional to the number of random 
steps in the time window (t — t* ) . In the same way, we 
may define x{t) for t > t* as 



x{t) = + (< - t*) J, 



(5) 



where again x* — x{t*) and J represents the steady (dc) 
current. Expressing {t — t*) from Eq. ([5]) and inserting it 
into Eq. (j4]), we may express time-dependence of d±{t) 
through x{t). To justify the choice of ansatz in Eq. (j4|), 
we introduce a quantity £^{t) defined as 



x{t) 



A 
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(6) 



where A = d\{t*) — ^x* (in our case, A < 0). For long 
enough times, we expect A/x{t) — > and (,{t) should 
approach the constant f ~ D±/j. Indeed, we observe 
in Fig. im^d) the saturation of ^(t) for different values 
of F and J. Moreover, for a fixed J the values of ^ are 
independent of the strength of electric field F. This result 
implies that for moderate driving the number of random 
steps {t — t*)D± is proportional to the distance travelled 
along the field {t — t*)j. This proportionality suggests 
yet another strong argument supporting the decay of the 
driven bound state in 2D. 

Decay of the bound state can be also monitored by 
calculating the 2D correlation function C(r) measuring 
time-dependent probability for the hole pair to be at a 
relative position r. Results in Fig. [3] show a disk-shaped 
pattern of C(r) elongated perpendicular to the field, con- 
sistent with Fig. [2][a). Moreover, a perpendicular cut 
through r = (0,0) (not shown) reveals that the position 
of maximum of C(r) steadily moves to larger |r|, deter- 
mining the main direction of decay. 



Discussion and Conclusion. — By applying the i-ED 
and <-EDLFS method to study real-time response of a 
fully quantum system, we managed to follow the out of 
equilibrium dynamics of a driven system where initially 
at i = the bound state exists due to the exchange of spin 
excitations. Our calculations on the ladder system show 
that as long as there is no additional mechanism to pro- 
vide the glue for binding, the bound state of two charge 
carriers decays with the onset of finite steady current. 
In the 2D system a bound pair decays predominantly in 
the perpendicular direction with respect to the external 
field, which consequently allows for more efficient release 
of the gained potential energy through magnon emission. 
At longer times, however, the motion of carriers perpen- 
dicular to the field can be consistently described by a 
random walk with the same scattering mechanism as for 
the propagation along the field. Therefore, assuming that 
preformed pairs exist in a superconductor above Tc [30| , 
our data indicate that steady propagation of bound pairs 
may not be realized as long as paring and dissipation em- 
anate from identical degrees of freedom. 

Our results on the decay of the bound state may as 
well lead to a broader understanding of driven strongly 
correlated systems. We found a significant difference in 
out-of-equilibrium response between quasi-lD and 2D 
systems. In this context, we showed that at short times 
after switching on the field the perpendicular distance 
between carriers dj_(t) is universally determined by the 
distance travelled by the center of mass of two carriers 
x{t). Due to additional decay channels that open as a 
consequence of charge motion along the perpendicular 
direction, the characteristic decay time on the square lat- 
tice is much shorter than on the ladder system. We ex- 
pect that such cooperative correlation between parallel 
and perpendicular dynamics may also manifest itself in 
various setups driven away from equilibrium (like that in 
Ref . [sH ) where charge carriers initially form a state with 
inhomogeneous microscopic structure. 
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